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ABSTRACT 

Due to the exquisite photometric precision, transiting exoplanet discoveries from the Kepler mission 
are enabling several new techniques of confirmation and characterization. One of these newly accessible 
techniques analyzes the phase variations of planets as they orbit their stars. The predicted phase 
variation for multi-planet systems can become rapidly complicated and depends upon the period, 
radius, and albedo distributions for planets in the system. Here we describe the confusion which may 
occur due to short-period terrestrial planets and/or non-transiting planets in a system, which can 
add high-frequency correlated noise or low-frequency trends to the data stream. We describe these 
sources of ambiguity with several examples, including that of our Solar System. We further show 
how decoupling of these signals may be achieved with application to the Kepler-20 and Kepler-33 
multi-planet systems. 

Subject headings: planetary systems - techniques: photometric 



1. INTRODUCTION 

The past few decades represent an extraordinary pe- 
riod of growth in our knowledge of planetary systems. 
Whereas previously we used the single data point of our 
own Solar System to derive theories of planetary forma- 
tion, we now have a diverse range of planetary systems to 
draw upon. We are gaining an improved understanding 
of the distribution of planetary properties such as pe- 
riod, mass, eccentricity, radius, and multiplicity. These 
last two have been aided in no small part by the signifi- 
cant discoveries of the Kepler mission which has revealed 
nume rous cases of transitin g planets in multi-pla net sys- 
tems (|Borucki et al.ll201 1a b: Bat alha et al1f2012h . Many 
of these are tightly-packed systems of multiple plan- 
ets within a relatively s mall period range, exa mples of 
whi ch include Kepler-11 dLissauer et alJl2011afl . Kepler- 
20 (iFressin et al.1120121; iGautier et al.ll2012f) and K epler- 
33 dLissauer et al.l I2012T) ■ iLissauer et al I (|2011b[ ) and 
iLissauer et al.l (|2012f) point out that these multi-planet 
systems comprise a large fraction of the total candidates 
detected and are far less likely to be due to false-positives. 

The precision of the Kepler photometry has also al- 
lowed the investigation of out-of-transit variations which 
are phased with the planetary orbit. These include ellip- 
soidal variations of the host star induc ed by the planet 
(jJackson et al.l l2012t IPfahl et al.ll2008D. Doppler boost - 
ing or beaming ([Mazeh et al. l l2012HShporer et al.llMTI). 
and reflecte d light from the planet ([Kane fc GelinoluOlOl . 
120111 I2012D . A combination of these effects have been 
detected in several cases, such as HAT-P-7b which was 
observed in Kepler data to have sig natures of reflecte d 
light as well as ellipsoidal variations ( Welsh et al.l f2010). 
These effects have almost exclusively been considered for 
single-planet systems, usually in very short period orbits 
where these effects will be higher in amplitude. 

Here we consider the confusion that can result from 
multi-planet systems when attempting to measure the 
photometric phase variations resulting from reflected 
light. Since there will be a variety of periods, radii, 
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and albedos, the contribution of each planet to the total 
phase variations may encompass a large range of am- 
plitudes depending on how these values are distributed. 
The result of this is that certain planetary configurations 
will have "rogue" phase variation components which 
masquerade as correlated noise in the data. We describe 
a decoupling procedure for these systems to disentangle 
the signatures of the individual planets. We discuss the 
implications for the Solar System as an externally ob- 
served exosystem and apply the method to the Kepler 
multi-planet systems of Kepler-20 and Kepler-33. 

2. PHOTOMETRIC PHASE VARIATIONS 

The phase variations of exoplanets are described in de- 
tail from numerous sources, such as those mentioned in 
Sec tion ffl We refer the re ader to iKane k Gelinol (|2010f l 
and lKane fc Gelinol ([20111 ) since we will be adopting that 
particular formalism. 

The flux ratio of an exoplanet to the host star as ob- 
served from Earth is given by the following expression 



e(a, A) 



" A(A) 



= A g (X)g(a,X)'- 



Rl 



(1) 



where A is the wavelength of the observations, A g (\) is 
the geometric albedo, g(a, A) is the phase function, and 
R p is the radius of the planet. The phase angle of the 
planet, a, is defined to be zero when the planet is at supe- 
rior conjunction. Here we also allow for orbital eccentric- 
ity by using the time-dependent star-planet separation, 
r, given by 

a(l-e 2 ) 



1 + e cos / 



(2) 



where a is the semi-major axis and e is the orbital ec- 
centricity. Notice that Equation [T] contains three major 
components: the geometric albedo, the phase function, 
and the inverse-square relation to the star-planet sep- 
aration. For a given phase function, the flux ratio is 
proportional to the albedo, the square of the planetary 
radius, and the inverse square of the separation. This is 
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Fig. 1. — The model photometric flux variations due to planetary 
in tightly-packed orbits close to the star and a Jupiter-size planet at 
planets and the dotted line indicates the combined effect. 

equivalent to being approximately proportional to P 3 / 4 
where P is the orbital period. In other words, 



e(a, A) oc A g x R 2 p x P 3 / 4 



(3) 



The radius is clearly the dominant component which in- 
fluences the amplitude of the phase variation. As we will 
see, this is an important factor when assessing the sig- 
natures of multi-planet systems, particularly the tightly- 
packed systems detected by the Kepler mission. 

3. MULTI-PLANET SYSTEMS 

Here we describe the combination effects and decou- 
pling issues for multi-planet systems. 

3.1. Combination Effects 

The discoveries from both RV and transit exoplanet de- 
tections have revealed that planetary systems can come 
in a variety of orbital configurations and physical char- 
acteristics. Multiplicity appears to be a common trait 
of exosystems and thus we can expect most systems to 
exhibit phase signatures of more than one planet. The 
manner in which these signatures contribute to the to- 
tal flux variations can cause an ambiguous interpretation 
of the data. Short-period planets can cause significant 
confusion if the observational cadence is insufficient to 
sample the variation cycle caused by the planets. For 
systems in which there is a dominant giant planet, the 
signatures of these smaller short-period planets can ap- 
pear as correlated noise in the data depending on the 
photometric precision. In addition, it is unlikely that all 
the planets in a given system will transit their host star 
from the perspective of Earth, as we will see for the So- 
lar System in Section [4.21 Thus there may be additional 
unaccounted for signatures in the data if one is basing 
the model on those known to transit. 

Shown in Figure [T] is an example model of flux varia- 
tions due to the phase signatures of a four-planet system, 
the contribution from each planet shown as a solid line. 
The dominant source of the reflected light is due to a 
Jupiter-size planet in a circular orbit at 0.7 AU. There 
are also three Earth-size planets with periods less than 
30 days. The combined signature is represented by a dot- 
ted line. Although one may detect the signal of the giant 
planet, the terrestrial planets will add correlated noise to 



phases in a hypothetical system consisting of three Earth-size planets 
0.7 AU. The solid lines show the phase variations due to the individual 



the signal, the amplitude of which will depend on their 
respective albedos. The opposite effect is also true, that 
the phase signature of three known (from their transits) 
terrestrial planets will soon be distorted by the presence 
of an unknown (non-transiting) exterior giant planet. 

3.2. Fourier Decoupling 

For multi-planet systems, the planets are more likely 
to be in near-c ircular orbits than planets in single- 
planet systems ([Wright et al.l [2009) . This circulariza- 
tion of the orbits is even more common for the tightly- 
packed planetary orbi ts found for the Kepler systems 
(|Moorhead et alll2011[ ). Thus, the combined phase sig- 
nature of multi-planet systems may be adequately de- 
scribed as a sum of trigonometric functions suitable 
for fourier analysis. Consideration of phase signa- 
tures for eccen t ric pla nets, such as those described by 
IKane k Gelinol ((2010), may also be considered in this 
context using Keplerian techniques in the fourier method. 
This is particularly important when evaluating the signif- 
icance of peaks in the powe r spectrum and is discus sed in 
detail bv ICummineJ (2004) . iCumming et al.l (|2008f> . and 
IQ'Toole et al.1 (120091) . 

The capability of a robust fourier analysis to recog- 
nize the individual component planetary contributions 
to the total phase variations depends primarily on (a) 
the photometric precision, (b) the observational cadence, 
and (c) the duration of the observations. The observa- 
tional cadence and duration limits the period analysis 
between the low fundamental frequency and the Nyquist 
frequency. A low signal-to-noise dataset will produce 
peaks in the fourier spectrum of similar power to spu- 
rious signals (aliases). An example is shown in Figure 
[2] where we have simulated a four-planet system with 
radii ranging from super-Earth (2R®) to large ice giants 
(5P©). The orbits are tightly packed with the outer 
planet in a ~ 25 day period orbit. The dotted line 
shows the combined phase variation for the entire sys- 
tem of planets. We have simulated data for 200 epochs 
spread over twice the period of the outer planet. Each 
time was passed through a gaussian filter with a stan- 
dard deviation of 15 minutes to randomize the times. 
A weighted Lomb- Scargle (L-S) period analysis ([Lombl 
Il976t IScarglelll982l ) was applied using the methodology 
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Fig. 2. — Top panel: The flux variations exhibited by a hypothetical closely- packed system of four super-Earths. The combined signature 
of all four planets is shown as a dotted line, along with simulated data acquired over the time-span of one complete orbit of the outer 
planet. Bottom panel: The periodogram resulting from a fourier analysis of the data. The vertical dashed lines indicate the actual periods 
of the four planets. 

ted signals and then re- analysing the residuals can re- 
veal the remaining planets. This method however does 
assume gaussian noise with purely periodic signatures. 
Since the number of cycles is very limited, the Maximum 
Entropy Method (MEM) is an alternative approach since 
this is relatively efficient in detecting frequency lines with 
few assumptions regarding t he initial estimate s of the fit 
parameters (see for example lEndl et all (|2000ft ). 



shown in lAharmin et all (|2005f ) . This method applies in- 
verse variance weighting calculated for each data point 
to each of the terms in the fourier component. In cases 
where the error bars (and hence the weights) are equal, 
the weighted L-S periodogram reduces to the standard 
L-S periodogram. The bottom panel shows the peri- 
odogram for the simulated data with vertical dashed lines 
indicating the actual periods of the four-planet system. 

There are a couple of issues to note regarding the 
Fourier disentanglement of the individual signatures de- 
picted here. The planet located at ~ 17 days has the 
lowest phase amplitude, comparable to the noise proper- 
ties of the data, and therefore is barely detectable. As 
such, improving either the cadence or the time baseline of 
the observations has little effect in improving the detec- 
tion of this signature. The detection of the planet with 
the smallest period is only affected in so far as the ca- 
dence drops to a level that the orbital frequency becomes 
comparable to the Nyquist frequency. In this case, that 
situation would result if the number of measurements 
dropped from 200 to 27 yielding a cadence of ~ 1.8 days. 
One must consider this detection threshold when plan- 
ning the observation strategy for a particular mission. 

The fourier method is thus suitable to extract such sig- 
natures in most cases except where the signal-to-noise is 
very low. Where one or more phase signatures dominate 
the signal, an iterative approach of subtracting those fit- 



4. THE SOLAR SYSTEM 

To find a classical case where phase variation confusion 
can occur, we need look no further than our own Solar 
System. The data for the Solar System planets were 
extracted from the JPL HORIZONS Systerrfl 

4.1. Combined Phase Variations 

First we consider the total phase variations expected 
by the Solar System when viewed along the ecliptic. Ta- 
ble Q] contains the Solar System orbital parameters used 
for these calculations. Note that when treating the Solar 
System as an exosystem, one must be careful to distin- 
guish between the longitude of perihelion and the ar- 
gument of perihelion since the longitude of perihelion 
is the sum of the longitude of the ascending node and 
the argument of perihelion. Here we use the argument 

1 http : //ssd . jpl . nasa . gov/?horizons 
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of perihelion (w) derived from these quantities. From 
these parameters and the geometric albedos, we calcu- 
late the predicted phase variations whose amplitudes are 
also shown in Table [TJ 

Figure [3] shows these phase variations at two resolu- 
tions; one phased on the orbit of Neptune and the other 
phased on the orbit of Jupiter. The top panel is domi- 
nated by Jupiter which has the largest phase amplitude. 
Uranus and Neptune have negligible contributions to the 
combined phase variations and Saturn produces a small 
modulation in the Jupiter signature. However, the ter- 
restrial planets are at least as significant in their con- 
tributions. Venus has a particularly large contribution, 
the amplitude of which is 43% that of Jupiter. The ef- 
fects of the terrestrial planets can be seen much clearer 
in the lower panel which shows the phase variations for 
one complete orbital period of Jupiter. Reaching the 
precision to detect the largest planet in our system is no 
guarantee that one would see the effects of the terrestrial 
planets as anything other than partially-correlated noise. 
This is compounded by the relatively short orbital peri- 
ods of the inner planets whose relation to the observation 
cadence may render those planets all but undetectable. 
The flux amplitudes for the Solar System planets are ex- 
tremely small, but this is an effect which can be scaled to 
alternative orbital configurations. We will see several ex- 
amples of this in the following section discussing Kepler 
systems, where the predicted flux amplitudes are several 
orders of magnitude larger. 

4.2. Transit Probabilities 

It is worth noting that not all of the Solar System 
would be detectable via the transit method, even if 
viewed along the ecliptic. The detection of the outer 
planets via either the transit or radial velocity techn iques 
is especially difficult, as described bv iKand (|2011f ). For 
a circular orbit the geometric transit probability is in- 
versely proportional to a, such that the inclination of 
the planet's orbital plane, i, must satisfy 

a cos i < R p + R± (4) 

By this criteria, observing our Sun along the ecliptic 
plane would only result in a transit of Earth. We per- 
formed a Monte-Carlo simulation which rotates the plane 
of the Solar System ±10° with respect to the observer. 
This simulation showed that there is only a very nar- 
row range of viewing angles for which one would see 
more than one planet transit the Sun. Specifically, if 
the viewing angle departs from the ecliptic by between 
1.76°and 1.78°, one will observe transits of both Mars 
and Neptune. This could be considered an extreme case 
of demonstrating that not all planets are necessarily ac- 
counted for in multi-planet transiting systems and thus 
the phase variations shown in Figure [3] will not be rep- 
resentative of what is predicted. However, the issue is 
still particularly pertinent for Kepler multi-planet sys- 
tems which are in much more compact orbital configura- 
tions, as we will see in the following sections. 

5. KEPLER SYSTEMS 

Here we apply this methodology to several of the 
known multi-planet systems detected by the Kepler mis- 
sion and discuss the impact on phase signatures. 



5.1. Kepler- 20 

The Kepler-20 system consists of five known transiting 
plan ets which wer e disco vered by iFressin et al.l ()2012l ) 
and iGautier et al.l (|2012l ). Three of the planets (b, c, 
and d) are of super-Earth size 1.91-3.07 i?® and the 
other two (e and f) are Earth size and smaller. The 
orbital and physical parameters for these planets are 
shown in Table [2j In the absence of information regard- 
ing the atmospheric properties of these planets, we have 
no knowledge of their geometric albedos. The poten- 
tial diversity of exoplanetary atmospheres, even amongst 
those of similar mass and semi-major axis, is sufficient 
to disqualify an assumption regarding the albedo of a 
planet. However, we are aware of an albedo depen- 
dence on star-planet separation dSudarskv et al.l 120051 : 
ICahov et al.lfe OlO: Ka ne fc Gelinoll2010D . such as the case 
of HD 209458b which was determined to have an up- 
per limit of A g < 0.08 from observations using the Mi- 
crovariability and Oscillations of STars (MOST) satellite 
(R owe et al]l2008l ). For the purposes of this analysis, we 
assign a moderately low default value of A g =0.2 to all 
planets so that we may evaluate the relative contribu- 
tions of the other parameters to the overall flux. When 
performing the analysis for a given dataset, these albe- 
dos may be treated as free parameters when identifying 
the shape of the phase signatures which match to the 
observed orbital periods from the transit photometry. 

The resulting model phase variations for the Kepler-20 
system are shown in the top panel of Figure 21 with the 
contribution from individual planets shown as solid lines 
and the total phase variations indicated by a dotted line. 
The two largest contributors to the phase variation are 
the b and c planets due to their larger size and prox- 
imity to the star, shown in Table [2] The e planet has 
a smaller star-planet separation than the c planet but 
the phase amplitude is inhibited by its small size. How- 
ever, were the e planet to have a substantial reflective 
atmosphere such as Venus, the amplitude of the phase 
variation for planet e would rise to become comparable 
to around half that of planet c. Thus the two smallest 
planets in this system may yet be distinguishable from 
the major sources of reflected light in the system. 

The bottom panel of Figure [4] shows a sample peri- 
odogram using the techniques described in Section 13.21 
where we have simulated data for the model phase curves 
shown in the top panel. The vertical dashed lines indi- 
cate the orbital periods of the known five planets in the 
system. As discussed earlier, in cases such as this where 
the phase signature is largely dominated by one or more 
planets, the remaining signatures may remain invisible 
regardless of the quality of the data or the time baseline 
of the observations. Thus we see that the signatures of 
the b and c planets, at 3.7 and 10.9 days, are unambigu- 
iously extracted from the data. One possible solution is 
to perform a fit to the these two phase signatures, sub- 
tracting them from the data, and then re-analysing the 
residuals in an attempt to reveal presence of the remain- 
ing planets. The success of such an approach depends 
on the amplitude of the remaining phase signatures with 
respect to the noise properties of the data. 

To investigate the potential for additional sources 
of phase amplitude components in the system, we re- 
analyse the Keck/HIRES radial velocity data acquired 
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TABLE 1 

Solar System Planetary Orbital Parameters and Peak Flux Ratios 



Planet 


P 


a 


e 


id 


R p 






Flux Ratio 




(days) 


(AU) 




(deg) 




(deg) 




(io- 9 ) 


Mercury 


87.97 


0.39 


0.206 


29.1 


0.38 


83.0 


0.11 


0.172 


Venus 


224.70 


0.72 


0.007 


55.2 


0.94 


86.6 


0.65 


2.000 


Earth 


365.26 


1.00 


0.017 


114.2 


1.00 


90.0 


0.37 


0.646 


Mars 


686.98 


1.52 


0.093 


286.5 


0.53 


88.2 


0.15 


0.040 


.Jupiter 


4332.82 


5.20 


0.048 


275.1 


11.22 


88.7 


0.52 


4.836 


Saturn 


10755.70 


9.54 


0.054 


336.0 


9.46 


87.5 


0.47 


0.880 


Uranus 


30687.15 


19.19 


0.047 


96.5 


4.01 


89.2 


0.51 


0.037 


Neptune 


60190.03 


30.07 


0.009 


265.6 


3.89 


88.2 


0.41 


0.013 




Orbital Phase 



Fig. 3. — Predicted photometric flux variation for the Solar System as observed along the ecliptic. The top panel is phased on the orbital 
period of Neptune. The bottom panel is phased on the orbit of Jupiter to highlight the contributions of the terrestrial planets, most notably 
Venus. The dashed lines in each case represent the normalized phase function of Neptune and Jupiter respectively. 



TABLE 2 

Kepler-20 Planetary Orbital Parameters and Peak 
Flux Ratios 



Planet 


P 


a 


R P 




A 3 


Flux Ratio 




(days) 


(AU) 


(Re) 


(deg) 




(10" 6 ) 


b 


3.696 


0.045 


1.91 


86.50 


0.2 


0.639 


c 


10.854 


0.093 


3.07 


88.39 


0.2 


0.394 


d 


77.612 


0.345 


2.75 


89.57 


0.2 


0.023 


e 


6.098 


0.063 


0.87 


87.50 


0.2 


0.068 


f 


19.577 


0.138 


1.03 


88.68 


0.2 


0.020 



for Kepler-20 bv lGautier et al.l ()2012l ). In particular, we 
would like to investigate evidence of any additional com- 
panions in the system which may not necessarily transit 



the host star. Shown in the left panel of Figure [3] are 
the RV data along with the Keplerian m odel predicted 
by the planetary parameters described bv lGautier et al.l 
(2012) and the residuals to the fit shown at the bottom. 

The mass of the d planet is difficu lt to constra i n from 
the RV data and is described by iGautier et al.l (|2012f ) 
as having a 2a limit of < 20.1 M®. The Markov Chain 
Monte Carlo (MCMC) analysis of the RV data yields a 
mass estimate for the d planet of 9.50lg'^ M® (Geoff 
Marcy, private communication). We adopt this value of 
the mass for deriving the Keplerian model for the data 
which then produces the shown residuals. The predicted 
RV semi-amplitudes for the e and f planets are 0.38 m/s 
and 0.47 m/s respectively and therefore have a negligible 
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Fig. 4. — Top panel: The model photometric flux variations due to planetary phases in the Kepler-20 system for one complete orbital 
period of the outer (d) planet. The solid lines show the phase variations due to the individual planets and the dotted line indicates the 
combined effect. Bottom panel: The periodogram resulting from a fourier analysis of the data. The vertical dashed lines indicate the actual 
periods of the five planets. 




Fig. 5. — Left panel: The radial velocity data of the Kepler-20 system (Gauticr ct al. 2012) with the radial velocity model predicted 
by the planetary parameters. The lower panel shows the residuals after subtracting the model from the data. Right panel: Regions of 
exclusion (above the lines) for additional planets within the system based upon the rms scatter of the residuals. The detected planets in 
the system are shown for reference. 
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TABLE 3 

Kepler-33 Planetary Orbital Parameters and Peak 
Flux Ratios 



Planet 


p 


a 


R P 


i 


A 3 


Flux Ratio 




(days) 


(AU) 


(«a) 


(dcg) 




(10" 6 ) 


b 


5.668 


0.068 


1.74 


86.39 


0.2 


0.238 


c 


13.176 


0.119 


3.20 


88.19 


0.2 


0.262 


d 


21.776 


0.166 


5.35 


88.71 


0.2 


0.375 


c 


31.784 


0.214 


4.02 


88.94 


0.2 


0.128 


f 


41.029 


0.254 


4.46 


89.17 


0.2 


0.112 



contribution to the model. 

The residuals shown in Figure [S] still contain a slightly- 
reduced rms scatter of 4.5 m/s. We calculate RV ampli- 
tudes that equal this la scatter as a function of orbital 
period and planet mass. These are shown in the right 
panel of Figure [5j with the lines representing multiples of 
the rms scatter of the residuals. One consequence of this 
figure is that we can exclude the presence of a Jupiter- 
mass planet in the system out to a period of ~ 1000 days 
at the 3cr level. The exclusion yields significant informa- 
tion for a correct model of the phase variations since a 
Jupiter-size planet with a 200 day period would produce 
a long-term trend in the photometry, where the ampli- 
tude is similar to that of the c planet. 

5.2. Kepler-33 

Kepler-33 is a system of five known transiting plan- 
ets which was announced bv lLissauer et all ((2012) . The 
complete parameters for these planets are shown in Table 
[21 where we have once again adopted fiducial albedo val- 
ues of 0.2. This system is quite different from the Kepler- 
20 system since it consists of planets in a size regime 
ranging from super-Earth to Neptune and a compact pe- 
riod range between 5 and 42 days. The major conse- 
quence of this compact planetary configuration is that all 
of the detected planets contribute comparable amounts 
of flux to the total phase variations of the system. For 
example, since the outer planets have a larger size, it 
somewhat compensates for their larger star-planet sepa- 
ration. The phase variations are shown in the top panel 
of Figure HI where the solid lines indicate the individual 
planets and the dotted line represents the combined ef- 
fect. The similar phase amplitudes and small range of 
periods results in a significantly improved fourier decou- 
pling of the phase signatures compared with the Kepler- 
20 system. The bottom panel of Figure [6] shows the peri- 
odogram from a simulated dataset constructed from the 
phase curves shown in the top panel. The inner three 
planets are recovered with ease due to their relatively 
large phase amplitudes. The outer two planets have de- 
tectable signatures from the combined analysis of all five 
planets but would benefit from the prior extraction of 
the three dominant planets, as discussed earlier. 

RV data is not available for Kepler-33 at the time of 
writing, which prevents the exclusion of additional plan- 
ets as was performed for the Kepler-20 system. For such 
a compact system, it is not unreasonable to postulate the 
existence of further non-transiting planets in the system. 
This could have a profound effect on the combined phase 
variations if Jovian-size planets lie just beyond the reach 
of the transit method to detect them. An example of such 
a case is shown in Figure [7] We have added the signature 



(shown as a dashed line) of a non-transiting Jupiter-size 
planet with an orbital period of 55 days. This planet 
immediately dominates the phase variations with an am- 
plitude which is almost a full order of magnitude larger 
than the previous largest amplitude of the c planet. The 
signatures of the inner planets are thus reduced to cor- 
related noise in the face of this new dominant signature, 
making extraction more difficult since the albedos of all 
of the planets are treated as free parameters in the fit to 
the data. 

6. CONCLUSIONS 

Initial results from the Kepler mission indicate that 
multi-planet systems are relatively common. The preci- 
sion of the Kepler mission and future photometric studies 
will allow detailed investigations of photometric varia- 
tions due to planetary phase signatures. We have shown 
that the various phase signatures in these systems can be 
combined to produce an emulation of correlated noise in 
the photometry. This confusion can both cause the plan- 
ets to be indistinguishable from each other and disguise 
the phase signature from the dominant planet in the sys- 
tem. On the other hand, non-transiting planets, usually 
at star-planet separations beyond the known size of the 
system, may either introduce new sources of correlated 
noise or indeed be the dominant source of variations in 
the system, such as was demonstrated for Kepler-33. 

We have shown that fourier analysis of the data will 
be able to extract many signatures from multi-planet sys- 
tems, but depends critically upon the cadence and dura- 
tion of the observations in order for the frequency distri- 
bution of the signatures to be properly sampled. In prac- 
tice, the largest hinderance will be that of the photomet- 
ric precision, which will prevent even high-frequency sig- 
nals from being detected. Our investigation of the Solar 
System demonstrates an example of confusion sources, 
where the signals of Jupiter and Venus may be extracted 
but the signatures of the remaining planets are unlikely 
to be interpreted as anything other than correlated noise. 
It is worth noting that the Solar System would not be de- 
tected as a multi-planet system due to the mis-alignment 
of the planetary orbital planes, with the exception of 
Mars and Neptune. We further demonstrate the difficul- 
ties in deciphering the phase variation for the Kepler-20 
and Kepler-33 systems. The interesting aspect of these 
systems, along with many other Kepler multi-planet sys- 
tems, is the compact nature of the orbits such that it is 
possible to observe complete orbital phases for all plan- 
ets within a reasonable time frame. One corollary of this 
is that there is great potential for a non-transiting giant 
planet to exist beyond the known planets in each system 
whose signature would significantly change the predicted 
phase variations of the system. 

We have focused solely on the photoemetric variations 
due to the phase functions of the planets. There are 
however numerous other effects, mentioned in Section [TJ 
which need to be accounted for before one can begin to 
examine the signals discussed here. With further ob- 
servations of the host stars and an extended baseline of 
precision photometry, it is hoped that we can eventually 
disentangle the various contributions to the observed flux 
from these most interesting systems. 
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Fig. 6. — Top panel: The model photometric flux variations due to planetary phases in the Kepler-33 system for one complete orbital 
period of the outer (f) planet. The solid lines show the phase variations due to the individual planets and the dotted line indicates the 
combined effect. Bottom panel: The periodogram resulting from a fourier analysis of the data. The vertical dashed lines indicate the actual 
periods of the five planets. 
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Fig. 7. — As for the top panel of Figure l6l except that an additional non-transiting planet has been added with a period of 55 days, 
indicated by the dashed line. The plot now shows one complete orbital period of this new planet. 
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